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INTRODUCTICW 


This  chapter  studies  the  algebra  and  calculus  of  vectors  and  matrices,  as 
specifically  applied  to  the  USAF  Test  Pilot  School  curriculum.  The  course  is 
a  prerequisite  for  courses  in  Equations  of  Motion,  r^amics.  Linear  Control 
Systems,  Flight  Control  Systems,  and  Inertial  Navigation  Systems.  The  course 
deals  only  with  applied  mathematics;  therefore,  the  theoretical  scope  of  the 
subject  is  limited. 

Ihe  text  begins  with  the  definition  of  determinants  as  a  prerequisite  to 
the  remainder  of  the  text.  Vector  analysis  follows  with  rigid  body  kinematics 
introduced  as  an  application.  The  last  section  deals  with  matrices. 

2.2  DETERMINANTS 

A  determinant  is  a  function  \ghich  associates  a  number  (real,  imaginary, 
or  vector)  with  every  square  array  (n  columns  and  n  rows)  of  numbers.  The 
determinant  is  denoted  by  vertical  bars  on  either  side  of  the  array  of 
numbers.  Thus,  if  A  is  an  (n  x  n)  array  of  numbers  v^ere  i  designates  rows 
and  j  designates  columns,  the  determinant  of  A  is  written 


2.2.1  First  Minors  and  Cofactors 

When  the  elements  of  the  i  **  row  and  column  are  removed  from  a  (n  x 
n)  square  array,  the  determinant  of  the  remaining  (n-1)  x  (n-1)  square  array 
is  called  a  first  minor  of  A  and  is  denoted  by  M.^.  It  is  also  called  the 
minor  of  a^  ^ .  The  signed  minor ,  with  the  sign  determined  hy  the  sum  of 
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the  row  and  column,  is  called  the  cofactor  of  a^^  and  is  denoted  by 

Exanple: 

^11  ®12  ®13 

“  hijl  *  ®21  ^23 


then, 


^2  2 

^13 

M 

M32  = 

^3  2 

^33 

^2  3 

Mix  = 


Also,  A^i  =  (-!)'■''  M,,  =  (+1)  Mil  and  =  (-l)^'"^  M32  =  (-1) 


2.2.2  Determinant  Expansion 

The  determinant  is  equal  to  the  sum  of  the  products  of  the  elements  of 
any  single  row  or  column  and  their  respective  cofactors;  i.e.. 


1a|  =  a-iA-i  +  a.^A.^  +  ...  +  a.„A.„ 


A.  =  E  a. .A.  . ,  for  any  single  i  row. 

in  .  -  1 3  1 D 


column. 


2. 2. 2.1  Expanding  a  2  x  2  Determinant.  Expanding  a  2  x  2  determinant  about 
the  first  row  is  the  easiest.  The  sign  of  the  cofactor  of  an  element  can  be 
determined  quickly  by  observing  that  the  sums  of  the  subscripts  alternate  from 
even  to  odd  vdien  advancing  across  rows  or  down  colimns,  meaning  the  signs 
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alternate  also.  For  exanple, 


the  signs  of  the  associated  cofactors  alternate  as  shovm, 


+ 


By  deleting  the  row  and  column  of  a^^,  we  find  its  cofactor  is  just  the 
element  a^j  [actually  (+1)  x  ^22^  a  2  x  2  array,  and  likewise  the  cofactor 

®i2  ^~®2i  ^  (—I)  X  ^21^’  two  products  of  the  two 

diagonals  gives  us  the  expansion  or  value  of  the  determinant. 


a. 


11 


All 


+  a. 


12 


Ai2  = 


11  ®2  2 


+  d. 


12 


(-321 )  “ 


^11  ^2  2 


-  a 


12  ^21 


This  sinple  exanple  has  been  shown  for  clarity.  Actual  calculation  of  a 
2x2  determinant  is  easy  if  we  just  remember  it  as  the  subtraction  of  the 
cross  multiplication  of  the  elements.  For  exanple. 


(+)  (-) 
3 

5 

2. 2. 2. 2  Expanding  a  3  x  3  Determinant. 


(8)(5)-(3)(6)  =  22 


A|  =  la,.  = 


''^11 

^21 

^31 


a 


12 


■^22 


■*3  2 


a. 


13 


2  3 


■*3  3 
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Expemding  |a|  about  the  first  row  gives 

|a1“  ®i2^2  ~ 


®22 

^2  3 

^21 

^23 

^21 

^2  2 

(+1) 

+ai2  (-1) 

+ai3  (+1) 

^3  2 

^3  3 

®31 

^3  3 

^31 

^3  2 

®11  ^®22  ®33  “  ®23  ^32^  ^2  ^^21  ®33  ®23  ®31^  ■•■  3  ^®21  ®3  2  “  ®2  2  ®3 1  ^ 


Expanding  eind  grouping  like  signs. 


1  ®2  2  ®3  3 


+  s. 


12 


®2  3  ®31  ®13  ®21  ®3  2 


^13  ^2  2  ®3 1 


^11  ®2  3  ®3  2 


12  ^21  ^3 

for  3x3 


-  a 

XX  4C  ^  ox 

close  inspection  of  the  last  equation  shows  another  method 
determinants  using  diagonal  multiplication.  If  the  first  two  columns  are 
appended  to  the  determinant,  six  sets  of  diagonals  are  used  to  find  the  six 
terms  above.  The  signs  are  determined  by  the  direction  of  the  diagonal 
shown  in  the  illustration. 


(  +  )  {  +  )  (  +  ) 


=  (2)(5)(3)  +  (-1)(4)(1)  +  (6)(3)(-2)  -  (6)(5)(1)  -  (2)(4){-2)  -  (-1)(3)(3) 
=  30  +  (-4)  +  (-36)  -  30  -  (-16)  ^  (-9)  =  -40  +  25  =  -15 
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The  quicker  methods  of  calculating  determinants  are  useful  for  the  two 
simple  cases  here.  The  row  expansion  method  will  be  more  usefvd  for 
calculating  vector  cross  products.  The  use  of  determinants  for  solving  sets 
equations  will  be  discussed  later  in  this  chapter  in  the  matrix 
section.  Determinants  will  also  be  used  in  solving  sets  of  linear 
differential  equations  in  Chapter  3,  Differential  Equations. 

While  the  general  tool  for  evaluating  determinants  hand  calculation  is 
sisple,  for  determinants  of  greater  size  the  calculations  are  lengthy.  A  5  x 
5  determinant  would  contain  120  terms  of  5  factors  each.  Evaluating  larger 
determinants  is  an  ideal  task  for  the  conputer,  and  standard  programs  are 
available  for  this  task. 

2.3  VECTOR  AND  SCftLftR  DEFINITICT^S 

In  general,  a  vector  can  be  defined  as  an  ordered  set  of  "n"  quantities 
such  as  <aj ,  a^ ,  ...,  ajj>.  in  TPS,  vector  analysis  will  be  limited  to 
two-and  three-dimensional  space.  Thus,  xT  +  yj  and  xT  +  yj  +  zk  are 
representations  of  vectors  in  each  space,  while  xi,  yj,  and  zk  are  referred 
to  as  conponents  of  the  vector. 

Physically,  a  vector  is  an  entity  such  as  force,  velocity,  or 
acceleration,  vhich  possesses  both  magnitude  and  direction.  This  is  the  iisual 
approach  in  applied  physics  and  engineering,  and  the  results  can  be  directly 
applied  to  courses  here  at  the  School. 

Almost  any  physical  discussion  will  involve,  in  addition  to  vectors, 
entities  such  as  volume,  mass,  and  work,  vdiich  possess  only  magnitude  and  are 
known  as  scalars.  To  distinguish  vectors  from  scalars,  a  vector  quantity  will 
be  indicated  by  putting  a  line  above  the  symbol;  thus  ¥,  v,  and  a  will  be  used 
to  represent  force,  velocity,  and  acceleration,  respectively. 

The  magnitude  of  vector  F  is  indicated  by  enclosing  the  symbol  for  the 
vector  between  absolute  value  bars,  iTj.  Graphically,  a  scalar  quantity  can 
be  adequately  represented  by  a  mark  on  a  fixed  scale.  To  represent  a  vector 
quantity  requires  a  directed  line  segment  vdiose  direction  is  the  game  as  the 
direction  of  the  vector  and  \diose  measured  length  is  equal  to  the  magnitude  of 
the  vector. 
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The  direction  of  a  vector  is  determined  by  a  single  angle  in 
dimensions  and  two  angles  in  three  dimensions,  angles  \dJose  cosines  are  called 
direction  cosines.  This  text  vd.ll  not  deal  directly  vdth  direction  cosines, 
so  no  example  is  necessary. 


2.3.1  Vector  Equality 

TWO  vectors  whose  magnitude  and  direction  are  equal  are  said  to  be  equal. 
If  two  vectors  have  the  same  length  but  the  opposite  direction,  either  is  the 
negative  of  the  other.  This  is  true  even  vAien  graphically  two  vectors  are  not 
{diysically  drawn  from  the  same  starting  point. 

A  vector  that  may  be  drawn  from  any  starting  point  is  called  a  free 
vector .  However,  v^en  applied  in  a  problem,  the  position  of  a  vector  may  be 
important.  For  instance,  in  Figure  2.1,  the  distance  of  the  line  of 
application  of  a  force  fran  the  center  of  gravity  of  a  rigid  body  is  critical 
if  calculating  moments,  although  the  actual  point  of  application  along  the 
line  isn't  critical. 


FIGURE  2.1.  MOMENT  CALCULATION 

For  other  applications,  the  point  of  action  as  well  as  the  line  of  action 
must  be  fixed.  Such  a  vector  is  usually  referred  to  as  bound.  The  velocity 
of  the  satellite  in  the  orbital  mechanics  problem  shown  in  Figure  2.2  is  an 
exanple  of  a  bound  vector. 


ORBITAL  PATH 


POINT  OF  ACTION 


LINE  OF  ACTION 


FIGURE  2.2.  EXAMPLE  OF  A  BOUND  VECTOR 
2.3.2  Vector  Addition 

Graphically,  the  sum  of  two  vectors  A  and  B  is  defined  by  the  familiar 
parallelogram  law;  i.e.,  if  A  and  B  are  drawn  from  the  same  point  or  origin, 
and  if  the  parallelogram  having  A  and  B  as  adjacent  sides  is  constructed,  then 
the  sum  A  +  B  can  be  defined  as  the  vector  represented  by  the  diagonal  of  this 
parallelogram  vAiich  passes  through  the  common  origin  of  A  and  Vectors  can 
also  be  added  by  drawing  them  "nose-to-tail . "  See  Figure  2.3. 


FIGURE  2.3..  AMBITION  OF  VECTORS 
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Graphically  from  Figure  2.3,  it  is  evident  that  vector  addition 
commutative  and  associative,  respectively. 


A  +  B  -  B  +  A  and  A  +  (B  +  C)  »  (A  +  B)  +  C 


2.3.3  Vector  Subtraction 

Vector  subtraction  is  defined  as  the  difference  of  two  vectors  A  and  B, 

A  -  i  »  A  +  (-i). 


\diere 


(-1)(B)  -  (-B) 


and  is  defined  as  a  vector  with  the  same  magnitude  but  opposite  direction. 
See  Figure  2.4.  This  introduces  the  necessity  for  vector-scalar  multipli¬ 
cation. 


FIGURE  2.4.  VECTOR  SUBTRACTION 


2.3.4  Vector-Scalar  Multiplication 

The  product  of  a  vector  and  a  scalar  follows  algebraic  rules.  The 
product  of  a  scalar  m  and  a  vector  A  is  the  vector  mA,  \diose  length  is  the 
product  of  the  absolute  value  of  m  and  magnitude  of  A,  and  vrtiose  direction  is 
the  same  as  the  direction  of  A,  if  m  is  positive,  and  opposite  to  it,  if  m  is 
negative . 


2.3.5  Unit  emd  Zero  Vectors 

Regardless  of  its  direction,  a  vector  vdiose  length  is  one  (unity)  is 
called  a  unit  vector.  If  a  is  a  vector  with  magnitude  other  than  zero,  then 
unit  vector  a  is  defined  as  a/la|,  vAiere  a  is  a  unit  vector  having  the  same 
direction  as  a  and  magnitude  of  one.  It  happens  that  the  conpsnents  of  a  unit 
vector  are  also  the  cosines  of  the  angles  necessary  to  define  the  direction. 
Unit  vectors  in  the  body  axis  coordinate  system  will  retain  the  bar  symbol; 
i.e. ,  i,  j  and  k. 

The  zero  vector  has  zero  magnitude  and  in  this  text  has  any  direction. 
It  is  notationally  correct  to  designate  the  zero  vector  with  a  bar, 

2.4  LRW5  OF  VECTOR  ~  SCALAR  ALGEBRA 

If  A,  B,  and  C  are  vectors  and  m  and  n  are  scalars,  then 


1. 

mA  =  Am 

Commutative  Multiplication 

2. 

m(nA)  = 

(mn)A 

Associative  Multiplication 

3. 

(m  +  n)A 

=  mA  + 

nA 

Distributive 

4. 

m(A  +  B) 

=  mA  + 

TCB 

Distributive 

These  laws  involve  multiplication  of  a  vector  by  one  or  more  scalars. 
Products  of  vectors  will  be  defined  later. 

These  laws,  along  with  the  vector  addition  laws  already  introduced, 
enable  vector  ec[uations  to  be  treated  the  same  way  as  ordinary  scalar 
algebraic  equations.  For  exanple, 

if  A  +B  =  C 


then  by  algebra 


A  =  C  -  B 
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2.4.1  Vectors  in  Coordinate  Systems 

The  right-handed  rectangular  coordinate  system  is  used  unless  otherwise 
stated.  Such  a  system  derives  its  name  from  the  fact  that  a  right  threaded 
screw  rotated  throu^  90**  in  the  direction  from  the  positive  x-axis  to  the 
positive  y-axis  will  advance  in  the  positive  z  direction,  as  shown  in  Figure 
2.5.  Practically,  curl  the  fingers  of  the  right  hand  in  a  direction  from  the 
positive  x-axis  to  the  positive  y-axis,  and  the  thumb  will  point  in  the 
positive  z-axis  direction. 


FIGURE  2.5.  RIOT-HANDED  COORDINATE  SYSTEM 


An  inportant  set  of  unit  vectors  are  those  having  the  directions  of  the 
positive  X,  y,  and  z  axes  of  a  three-dimensional  rectangular  coordinate  system 
and  are  denoted  i,  j,  and  k,  respectively,  as  shown  in  Figure  2.5. 

Any  vector  in  three  dimensions  can  be  represented  with  initial  point  at 
the  origin  of  a  rectangular  coordinate  system  as  shown  in  Figure  2.6.  The 
perpendicular  projection  of  the  vector  on  the  axes  gives  the  vector's 
conponents  on  the  axes.  Multiplying  the  scalar  magnitude  of  the  projection  by 
the  appropriate  unit  vector  in  the  direction  of  the  axis  gives  a  conponent 
vector  of  the  original  vector.  Note  that  summing  the  conponent  vectors 
graphically  gives  the  original  vector  as  a  resultant. 


2.10 


FIGURE  2.6.  COMPONENTS  OF  A  VECTOR 

In  Figure  2.6  the  ccmponent  vectors  are  A^^i,  A^j,  and  Ajk.  The  sum  or 
resxiltant  of  the  ccm^nents  gives  a  new  notation  for  a  vector  in  terms  of  its 
conqponents. 


A  =  A^i+A2j+A3k 

After  noticing  that  the  coordinates  of  the  end-point  of  a  vector  A  vhose  tail 
is  at  the  origin  are  equal  to  the  components  of  the  vector  itself  (A^  =  x,  A^ 
»  y,  and  Aj  »  z),  the  vector  may  be  more  easily  written  as 

A  =  xi  +  yj  +  zk 

Hie  vector  from  the  origin  to  a  point  in  a  coordinate  system  is  called  a 
position  vector,  so  the  vector  notation  edaove  is  also  the  position  vector  for 
the  point  P.  The  same  definitions  for  notation,  components,  and  position  hold 
for  a  two-dimensional  system  with  the  third  component  eliminated. 

The  magnitude  is  easily  calcxilated  as. 


-Vv 


or 


2.11 


An  arbitrary  vector  from  initial  point  P(Xj 
vectors,  also. 


Yj ,  Zj )  and  terminal  poiri 


Zj )  such  as  shown  in  Figure  2.7  can  be  represented  in  terms  of  wit 


FIGURE  2.7.  ARBITRARY  VECTOR  REPRESENTATIC»I 


First  write  the  position  vectors  for  the  two  points  P  and  Q. 


and 


Cl  -  Xji  +  yj  +  Zjk 


Tj  =  x^i  +  y^  j  +  z^k 


Then  using  addition. 


or 


r^  +  PQ  = 

pQ  “=  >^2  “  +  (y2  “  yi  )T  +  (Zj  “  ^1 
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2.4.2  Dot  Product 

In  addition  to  the  product  of  a  scalar  and  a  vector,  two  other  types  of 
products  are  defined  in  vector  analysis.  The  first  of  these  is  the  dot,  or 
scedar  product,  denoted  by  a  dot  between  the  two  vectors.  The  dot  product  is 
an  operation  between  two  vectors,  and  results  in  a  scalar  (thus  the  name 
scalar  product).  Analytically,  it  is  calculated  by  adding  the  products  of 
like  ccmpjnents.  This  is,  if 


and 


then 


A  =  a^i  +  aj  j  +  ajk 


B  =  bji  +  bj  j  +  bjk 


A«B  =  a^bj  +  ajbj  +  ajbj 


vhich  is  a  real  number  or  scalar. 

Geometrically,  it  is  equal  to  the  product  of  the  magnitudes  of  two 
vectors  and  the  cosine  of  the  angle  between  them  (the  angle  is  measured  in  the 
plane  formed  by  the  two  vectors,  if  they  had  the  same  origin).  The  dot 
product  is  written 

A  •  B  =  |a1  jijcos  0 

Example 

{2i  -  ij  +  4k)  •  (-1  +  3 j  +  5k)  =  (2)(-l)  +  (-1)(3)  +  (4)(5)  =  15 

The  magnitudes  are 


and 


-\/4  +  1  +  16  =  4.6 
■V  1  +  9  +  25  =  5.9 


Therefore,  solving  for 

cos  0  =  15  /  (4.6)  (5.9)  =  15/27.1  =  0.553 

so 

0  =  56.1“ 
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Some  interesting  applications  of  the  dot  product  are  the  gecxnetrll 
implications.  For  instance,  the  geometric,  scalar  projection  of  one  vector  cm 
another  is  shown  on  Figure  2.8. 


lAlcose 

FIGURE  2.8.  GEOMETRIC  PROJECTION  OF  VECTORS 

Using  trigonometry,  the  projection  of  A  on  B  is  seen  to  be  equal  to  |a1  cos  6. 

A  quick  method  to  calculate  such  a  projection  without  knowing  the  angle  is  to 
cedculate  the  dot  product  and  divide  the  magnitude  of  the  vector  projected^ 
on.,  to.  •  That  is,  the  projection  of  A  on  to  B  is  equal  to  A  •  B/|b|  «=  |a1cos 

Several  particular  dot  products  are  worth  mentioning.  If  one  of  tl^^ 
vectors  is  a  unit  vector,  the  dot  product  becomes 


i  •  B  ■*  |il  1b|  cos  ©  =  (1)  |bI  cos  ©  «  IbI  cos  ©, 

^ich  is  the  projection  of  B  on  i  or  more  inportantly  the  conponent  of  B  in 
the  direction  of  i.  Also  note  the  dot  product  of  a  vector  with  itself  is  just 
equal  to  the  magnitude  squared,  since  the  angle  is  zero  and  cos  ©  *=  1.  More 
useful  is  the  situation  ^diere  two  non-zero  vectors  are  perpendicular 
(orthogonal).  The  dot  product  is  zero  because  the  cosine  of  90  degrees  is 
zero.  Thus,  for  non-zero  vectors  the  dot  product  may  be  a  test  of 
orthogonality.  Exanples  of  these  properties  using  standard  xmit  vectors  are 


1*1 


1  •  3 


j  •  k 
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2.4.3  Dot  Product  Laws 

If  A,  B,  and  C  are  vectors  and  m  is  scalar,  then 


1. 

A  •  B  «  B  •  A 

Commutative  Product 

2. 

A*{B  +  C)  >=A»B  +  A*'C 

Distributive  Product 

3. 

m(A  •  B)  =  (mA)  •  B  =  A  •  (itB) 

Associative  Product 

2.4.4  Cross  Product 

The  third  type  of  product  involving  vector  operations  is  the  cross,  or 
vector  product,  denoted  by  placing  an  "X"  between  two  vectors.  By  definition 
the  cross  product  is  an  operation  between  two  vectors  vdiich  results  in  another 
vector  (thus,  vector  product).  Again  both  analytic  and  geometric  definitions 
are  given. 

^^^lyhically,  the  cross  product  is  calculated  for  three-dimensional 
vectors  by  a  top  row  expansion  of  a  determinant. 


A  X  B  = 


a. 


+(-l) 


^3 

bs 

j  + 

(a^  bj  -  aj  bj  )i  +  (aj  b^  -  a^  )  j  +  (a^  b^  -  a^  b,  )k 


For  exanple,  _  _  _  _  _  _ 

(2i  +  4j  +  5k)  X  (3i  +  j  +  6k)  = 

1  I  k 

2  4  5 

3  16 

=  [(4)(6)-(5)(l)]i  -  [(2)(6)-(3)(5)]j  +  [ (2) (l)-(3) (4) jk 

=  19i  +  31  -  loic 
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The  geometrical  definition  has  to  be  approached  carefully  because  it  inusl 
be  remenbered  that  the  geometrical  definition  is  not  a  vector.  The  magnitude 
(a  scalar)  of  the  cross  product  is  equal  to  the  product  of  the  two  magnitudes 
and  the  sine  of  the  angle  between  the  two  vectors.  Thus 

|a  X  b|  •  |a|  |i|  sin  e 

While  the  magnit\ide  is  determined  as  above,  the  direction  of  the  resultant 
cross  product  vector  is  always  orthogonal  to  the  plane  of  the  crossed  vectors. 
Ihe  sense  is  such  that  vhen  the  fingers  of  the  right  hand  are  curled  from  the 
first  vector  to  the  second,  through  the  smaller  of  the  angles  between  the 
vectors,  the  thumb  points  in  the  direction  of  the  cross  product  as  shown  in 
Figure  2.9.  Note  the  importance  in  the  order  of  writing  A  X  B  since 
A  X  B  ^  B  X  A.  That  A  X  B  «  -B  X  A  is  easily  seen  using  the  right-hand  rule. 


FIGURE  2.9.  GEC»!ETRIC  DEFINITION  OF  THE  (3R0SS  PRODUCT 
The  cross  product  vector  u  can  be  represented  as 
U  *  A  X  B  *  I A I  I B 1  sin  0  u 

\here  u  is  a  unit  vector  in  the  direction  of  U,  ^ich  is  perpendicular  to  the 
plane  containing  A  and  B. 

Some  practical  applications  of  the  above  definitions  using  the  sine  of 
zero  and  90®  are  shown  for  \init  vectors  of  a  rectangular  coordinate  system. 


iXi  =  jXj  =  kxk  =  0  (note  the  zero  vector  has  any  direction) 

i  X  j  =  k  and  j  X  k  =  i  and  k  X  i  =  j  and 

j  X  i  =  -k  and  k  X  j  =  -i  and  i  X  k  =  -j 


These  cross  products  are  used  often,  and  an  easy  way  to  remeinber  them  is  to 
use  the  aid 


vdiere  the  cross  product  in  the  positive  direction  from  i  to  j  gives  a  positive 
k,  and  to  reverse  the  direction  gives  a  negative  answer, 

2.4.5  Cross  Product  Laws 

If  A,  B,  and  C  are  vectors  and  m  is  a  scalar,  then 

Anti-Commutative  Product 
Distributive  Product 
Associative  Product 

2.4.6  Vector  Differentiation 

The  following  treatment  of  vector  differentiation  has  notation  consistent 
with  later  courses  and  has  been  highly  specialized  for  the  USAF  Test  Pilot 
School  curriculum.  The  scalar  definition  of  the  time  derivative  of  a  scalar 
function  of  the  variable  t  is  defined  as, 

df(t)  _  lim  [  f(t  +  at)  -  f(t)  ] 
dt  At-»0  At 


1.  A  X  B  =  -B  X  A 

2.  ax(¥+c)=ax1b  +  axc 

3.  m(A  X  B)  =  (mA)  X  B  =  A  X  (mB) 
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Before  proceeding,  a  vector  function  is  defined  as 


F(t)  »  f^(t)i  +  fy(t)j  +  f^(t)k, 

\diere  f  ,  f  ,  and  f  are  scalar  functions  of  time  and  i,  j,  and  k  are  unit 
vectors  parallel  to  the  x,  y,  and  z  axes,  respectively.  A  vector  function  is 
a  vector  that  changes  magnitude  and  direction  as  a  function  of  time  and  is 
referred  to  as  a  position  vector .  It  gives  the  position  of  a  particle  in 
space  at  time  t.  The  trace  of  the  end  points  of  the  position  vector  gives  the 
trajectory  of  the  particle.  The  time  derivative  of  a  vector  function  with 
respect  to  some  reference  frame  is  defined  as. 


dF(t)  _  lim  [  F(t  +  At)  -  F(t)  ]  ^ 
dt"  ~  At^  At 


df  (t)  _ 
^  • 

“at —  ^ 


df  (t)  ^ 
-at—  3  + 


dfjt) 

“at 


df.  - 


df,  - 


df  _ 


“  ar  ^  ~3r  ^  ar 


k  = 


f,i  +  fyj  +  f,k 

vdiere  the  lack  of  a  function  variable  indicates  the  function  has  the  same 
variable  as  the  differentiation  variable,  and  the  dot  denotes  time 
differentiation. 


2.4.7  Vector  Differentiation  Laws 

For  vector  functions  A(t)  and  B(t),  and  scalar  function  f(t) 


1. 

d(A  +  B)  _ 
dt 

dA  dB 

at  dt 

Distributive  Derivative 

2. 

d(A  •  B)  _ 
dt 

—  dB  dA  -= 

^  •  at  +  at  •  ® 

Dot  Product  Derivative 

3. 

d(A  X  B) 

— at~ 

-  dB  dA  - 

^^at  ■'at*® 

Cross  Product  Derivative 

4. 

^  [f(t)B] 

dB  .  df(t)  -g 

-  at  "a^  ® 

Scalar,  vector  Product 
Derivative 
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2.5  LINEAR  VELOCIIY  AND  ACCELERATION 


Hie  time  derivative  of  a  position  vector  relative  to  some  reference 
system  is  the  linear  velocity.  Note  in  particular  that  the  velocity  of  a 
particle  is  a  vector  that  has  a  direction  and  a  magnitude.  The  magnitude  of 
the  velocity  is  referred  to  as  speed.  The  second  derivative  is  the  linear 
acceleration. 

Graphically,  the  derivative  of  a  vector  is  illustrated  as  shown  in  Figure 

2.10. 


FIGURE  2.10.  ILLUSTRATIW  OF  THE  DERIVATIVE  OF  A  POSITION  VECTOR 

The  difference  between  position  vectors  r(t  +  At)  and  r(t)  is  the 
numerator  of  the  definition  of  the  derivative.  The  arc  length  of  the 
trajectory  for  some  At  is  As.  If  we  neglect  the  division  At  and  are 
concerned  only  with  direction  of  the  derivative,  the  difference  of  the  two 
vectors  is  just  Ar  \frfiich  would  have  the  direction  as  shown  in  Figure  2.10. 
The  derivative  for  a  vector  r(t)  can  be  expanded  by  multiplying  by  the 
quantity  As/As  »  1,  as  follows. 


dr  ^  lim  ^  _  lim  ^  ^  lim  Ar  As 

Ht  “  At-K)  At  “  At-K)  At  As  “  At^  As  At 

but  as  At-K),  |Arl*  As,  therefore  lim  Ar/As  «  e  ,  since  its  magnitude  is  one 

At-K) 
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and  its  direction  becomes  tangential  to  the  trajectory.  The  As/At  portion 
gives  the  magnitvide  of  the  derivative  which  should  be  noted  as  the  speed  of 
the  particle  or  a  change  of  distance  per  time.  In  sxjmmary,  the  first 
derivative  of  a  vector  function  is  tangential  to  the  trajectory  and  has  a 
magnitude  that  is  the  speed  of  the  particle. 

Using  differentiation  law  four  to  take  the  derivative  of  the  vector 
written  in  the  form  of  magnitude  times  a  unit  vector, 

7(t)  =  r(t)r,  as  follows, 

dr(t)  d[r(t)r]  _  dr(t)  ...  dr 

"cE"  '  — at —  -  “at-  at 

note  that  the  linear  velocity  using  this  form  of  a  vector  has  two  components, 
the  first  is  the  rate  of  change  of  the  scalar  function  with  direction  the  same 
as  the  original  vector  itself.  The  second  conponent  is  the  scalar  function 
itself  with  the  rate  of  change  of  the  mit  vector  as  its  direction.  We  know 
that  the  roiit  vector  doesn't  change  magnitude,  but  it  may  change  direction 
giving  a  non-zero  derivative.  In  the  development  of  the  derivative  earlier, 
this  was  overlooked  since  the  rate  of  change  of  the  i,  j,  and  k  vectors 
are  fixed  in  a  coordinate  system  do  not  change  direction  or  magnitude. 

2.6  REFERENCE  SYSTEMS 


Linear  velocity  and  acceleration  have  meaning  only  if  expressed  (or 
implied)  in  reference  to  another  point  and  only  if  relative  to  a  particular 
frame  of  reference.  In  this  text  for  discussions  of  single  reference  systems, 
the  linear  velocity  and  acceleration  will  always  be  relative  to  the  origin  of 
the  reference  frame  in  vAiich  the  problem  is  given  and  will  be  denoted  by 
single  letters,  V  and  "i.  If  there  are  two  reference  systems  in  the  problem, 
the  notation  will  be  changed  to  read 


\diich  means  the  velocity  of  point  or  reference  A  relative  to  reference  B.  To 
take  a  time  derivative  of  a  vector  relative  to  reference  system  "A, "  the 
,  notation  will  be 
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^di 

“at 

There  shorild  be  less  confusion  in  multiple  reference  system  problems  con¬ 
cerning  \^ich  reference  frame  the  derivative  is  taken  by  using  this  notation. 
By  introducing  the  concept  of  multiple  reference  systems,  it  is  appropriate  to 
discuss  the  chain  rule.  For  two  reference  systems,  the  chain  rule  is  sinply 
stated.  For  point  A  in  reference  system  B,  ^diich  in  turn  is  in  another 
reference  system  C,  the  velocity  of  A  relative  to  C  is  equal  to 


A/C 


=  V, 


A/B 


+  V. 


B/C 


(2.1) 


While  calculating  derivatives  \dien  given  the  time  function  of  the 
trajectory  is  seemingly  single,  at  times  the  derivatives  may  be  difficult. 
Also,  if  the  fmction  is  not  known,  the  measurements  available  to  determine 
the  trajectory  may  be  in  terms  of  translational  or  rotational  parameters  idiich 
don't  always  lend  themselves  directly  to  a  time  function.  Another  method  of 
determining  velocities  and  accelerations  will  be  determined  using  pure 
translation  and  rotation.  Sinplification  will  consist  of  very  specific 
problems  with  convenient  alignment  of  reference  systems  at  specific  instances 
in  time.  So  it  will  appear  that  the  time  element  has  disappeared  in  the 
following  analysis  since  the  vectors  will  be  constants  at  the  instant  we 
observe  them. 


2.7  DIFFERENTIATION  OF  A  VECTOR  IN  A  RIGID  BODY 

Hie  two  basic  motions,  translation  and  rotation,  will  be  applied  to  a 
rigid  body  vdiich  is  assumed  not  to  bend  or  twist  (every  point  in  the  body 
remains  an  equidistance  from  all  others).  It  will  become  iitportant  to 
determine  not  only  the  velocity  and  acceleration  of  a  point  in  a  rigid  body, 
but  also  that  of  a  vector  vdiich  lies  in  a  rigid  body. 

2.7.1  Translation 

If  a  body  moves  so  that  all  the  particles  have  the  same  velocity 
relative  to  some  reference  at  any  instant  of  time,  the  bocty  is  said  to  be 
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in  pure  translation.  A  vector  in  pure  translation  cheuiges  neither  i 
magnitude  nor  direction  ^ile  translating,  so  its  first  derivative  would  be 
zero.  An  exanple  would  be  a  vector  from  the  center  of  gravity  to  the  wingtip 
of  an  airplane  in  straight  and  level,  unaccelerated  flight  with  respect  to  a 
reference  system  attached  to  the  earth's  surface.  From  the  ground  it  changes 
neither  magnitude  nor  direction,  although  every  point  on  the  aircraft  is 
traveling  at  the  same  velocity.  See  Figure  2.11. 


FIGURE  2.11.  TRANSLATION  AND  ROIATION  OF  VECTORS  IN  RIGID  BODIES 


2.7.2  Rotation 

If  a  body  moves  so  that  the  particles  along  some  line  in  the  body  have  a 
zero  velocity  relative  to  some  reference,  the  body  is  said  to  be  in  pure 
rotation  relative  to  this  reference.  The  line  of  stationary  particles  shown 
in  Figure  2.11  is  called  the  axis  of  rotation.  A  free  vector  that  describes 
the  rotation  is  called  the  angular  velocity,  w,  and  has  direction  determined 
by  the  axis  of  rotation,  using  the  right-hand  rule  to  determine  the  sense. 
The  chain  rule  as  described  for  linear  velocity  applies  to  the  angular 
velocity,  as  does  a  definition  of  its  magnitude  being  angular  speed.  The 
first  derivative  of  the  angular  velocity  is  the  angular  acceleration. 

It  can  be  proven  that  the  linear  velocity  V  of  any  point  in  a  rigid  body 
described  by  position  vector  r  \diose  origin  is  along  the  axis  of  rotation  can 
be  written 


2.22 


r 


V 


(0  X  r 


(2.2) 


Note  the  conventicms  using  the  right-hand  rule  apply,  and  V  is  perpendicular 
to  the  plane  of  T  and  «o. 

The  pure  rotation  of  one  reference  system  with  respect  to  another  would 
require  a  transformation  of  unit  vectors  frcan  one  system  to  another,  unless 
the  reference  systems  were  conveniently  aligned  at  the  instant  in  question. 
Such  transformations  are  considered  beyond  the  scope  of  this  course. 

Equation  2.2  can  be  generalized  to  include  any  vector  in  a  rigid  bo^ 
with  pure  rotation.  Refer  to  Figure  2.12. 


RIGID 

BODY 


FIGURE  2.12.  DIFFERENTIATION  OF  A  FIXED  VECTOR 

Let  p  be  a  vector  fixed  euiyvdiere  in  the  rotating  rigid  body  shown  in 
Figure  2.12.  The  problem  is  to  find  the  time  rate  of  change  of  the  vector. 
Two  position  vectors,  r^  and  r^ ,  from  the  origin  to  the  end  points  of  the 
vector  ^  are  drawn.  From  vector  addition 

■?!  +  p  =  I2 
P  «  r,  -  ii 

•  e  » 

P  '  -  r, 


or  solving 


Differentiating 
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From  Equation  2.2, 


Therefore,  the  derivative  of  any  fixed  vector  in  a  purely  rotating  rigid  body 
is  represented  the  cross  product  of  the  angular  velocity  of  the  rotating 
body  and  the  fixed  vector. 

2.7.3  Combination  of  Translation  and  Rotation  in  One  Reference  System  a 

It  is  possible  to  combine  the  two  types  of  velocity.  An  inportant  poin" 
to  notice  here  is  that  the  velocities  and  accelerations  are  arrived  at 
directly  without  the  use  of  position  vectors. 
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The  velocity  of  point  a  in  reference  system  D,  Figure  2.13,  will  be 
calculated.  The  rigid  bocfy  has  a  pure  angular  velocity,  w,  and  a  pure 
translation,  v,  in  reference  frame  D.  The  requested  velocity  is  just  the  sum. 


V  =  V  +  V 

rotation  translation 


Vrot. tion  is  equal  to  »  X  p.  is  given  as  v,  so 

V  =  «  X  p  +  V 

When  working  in  one  reference  system,  the  acceleration  may  be  calculated 
by  taking  the  derivative  of  the  velocity. 


TT  _  dv  _  d((o  X  p)  dv 
A  =  + 


~3t 


^  =  wXp  +  wXp 


+  V 


Here,  the  p  is  equal  to  to  X  p  as  was  shown  in  Ec[uation  2.3  and  "v  is  the 
translational  acceleration  a.  The  angular  acceleration  will  not  receive  any 
special  notation  in  this  text. 

So,  the  acceleration  in  a  single  reference  system  can  be  written 


A  =  wX(<oXp)  +  wXp  +  a 


2.7.4  Vector  Derivatives  in  Different  Reference  Systems 

The  more  general  problem  of  relative  motion  between  a  point  and  a 
reference  system  that  is  itself  moving  relative  to  another  reference  system 
will  be  approached.  More  than  one  reference  system  is  often  used  in  order  to 
sinplify  the  analysis  of  general  problems.  As  a  first  step,  it  is  necessary 
to  examine  the  procedure  of  differentiation  with  respect  to  time  in  the 
presence  of  two  references  moving  relative  to  each  other. 

A  reference  system  is  a  non-deformable  system  and  may  be  considered  a 
rigid  body.  So,  the  work  done  so  far  applies  here.  Figure  2.14  gives  the 
vectors  used  in  the  following  analysis. 
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FIGURE  2.14.  MOTION  WITH  TWO  REFERENCE  SYSTEMS 


Bie  problem  cdxnre  shows  point  p  with  position  vector  r  ,  moving  with 


P/B' 

respect  to  the  reference  B,  and  the  origin  of  B  with  position  vector  r 


B/C  ' 


moving  with  respect  to  reference  system  C.  The  reference  system  B  also  has  an 
angular  velocity  with  respect  to  C  of  The  goal  of  the  following 
development  will  be  to  find  the  time  rate  of  change  of  the  position  vector  in 
the  B  frame  as  seen  from  the  C  frame  or  notationally 


3t  *^P/B 

It  is  very  inportant  to  realize  that  this  is  not  the  same  as  the  velocity  of 
the  point  as  seen  from  the  C  frame.  Rather  it  is  the  rate  of  change  of  a 
position  vector  in  one  frame  as  seen  frcmi  another  frame.  So  the  derivative 
sought  is  not  This  velocity  would  be  obtained  using  the  chain  rule 
as  given  in  Equation  2.1. 

A  representative  exanple  is  the  motion  of  a  point  on  an  aircraft  with  a 
body  axis  system  at  the  center  of  gravity  and  the  aircraft  moving  along  some 
path  relative  to  the  ground.  The  second  reference  system  is  attached  to  the 
ground.  It  will  be  assumed  in  this  analysis  that  the  two  reference  systems 
have  the  same  unit  vectors.  Careful  attention  will  be  given  to  circumstances 
resrdlting  frcmi  axes  that  may  not  be  conveniently  aligned  during  the  analysis. 


Beginning  with  the  position  vector  in  the  frame  B, 


h/B  =  xi  +  yj  +  zk 

Differentiating  this  vector  with  respect  to  time  relative  to  the  C  reference 
frame  presents  a  problem  since  the  imit  vectors  of  the  B  system  are  rotating 
as  seen  in  the  C  system. 

So,  the  derivative  must  be  done  in  two  parts  using  the  fourth  law  given 
earlier. 


a?  S/B  “  xi  +  yj  +  zk  +  xi  +  yj  +  zk 

But  the  unit  vector's  derivatives  Ccin  be  written  as  vector  in  a  single 
reference  system  with  derivatives  as  seen  in  Equation  2.3.  Thus, 

i  =  X  i,  etc. 


So 

i  i 

at  S/B  *  xi  +  yj  +  zk  +  xi  +  yj  +  zk 

=  xi  +  ^  +  zic  +  x(<«^^^.  X  i)  +  y(oc^/c  X  7)  +  X  k) 

=  xi  +  ^  +  zic  +  X  (xi)  +  X  (yJ)  +  X  (zk) 

=  xi  +  ^  +  zic  +  X  (xi  +  yj  +  zk) 

The  first  three  terms  are  recognized  as  the  velocity  of  p  in  the  B  system  and 
the  next  term  is  the  cross  product  of  the  angular  velocity  of  the  B  system 
with  respect  to  the  C  system  and  the  position  vector  in  the  B  system. 


^  r 

Ht  P/B 


V.  * 


,  X  r  ,, 

C  p/B 


^p/B 


y,-  X  r 

/C  p/B 


(2.4) 
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This  equation  may  be  generalized  to  einy  vector  in  one  reference  systenf 
relative  to  another.  This  is  a  very  inqportant  relationship  and  will  be  used 
in  Chapter  A,  in  the  derivation  of  the  aircraft  equations  of  motion. 

The  acceleration  of  a  particle  at  point  p  would  be  handled  using  the 
definition  of  acceleration. 


p/c 


Ch  _ 

■r  V 
3t  ^P/c 


(2.5) 


Note  Equation  2.5  does  not  address  _ 

at  '^p/B 

Hopefully,  the  velocity  would  be  written  in  a  single  form  allowing  simple 
differentiation  to  obtain  the  acceleration.  If  not,  a  single  exchange  of 
notation  with  Equation  2.4  would  be  necessary. 


CX»!MENT  The  material  presented  thus  far  is  sufficient  to  enable  solution  of 
any  linear  or  angular  velocity  or  acceleration  in  a  kinematics 
problem.  However,  another  analysis  follows  which  may  clari 
multi-reference  problems  and  will  provide  definition  of  soi 
terms  that  will  be  of  value  in  later  courses. 


2. 7. 4.1  Transport  Velocity.  In  this  analysis  of  motion  relative  to  two 
reference  systems,  a  different  approach  is  taken  to  the  problem.  Figure  2.14 
is  expanded  as  shown  in  Figure  2.15  to  include  the  position  vector  directly 
from  reference  C  to  the  point  p. 


FIGURE  2.15.  TWO  REFERENCE  SYSTEM  VECTORS 


Ihus 


P/c 


’■p/B  ^B/C 


emd 


p/c 


*^p/c  at  S/C 


'-d  -  ''d  - 

dt  S/B  Ht  S/c 


\diere  the  first  term  is  Equation  2.4  and  the  second  is  V, 
these  terms,  _  _  _  _ 


p/c 


p/c 


S/B  +  ‘^/C  X  Ip/B  +  S/C 


.  S\]bstituting 
“  (2.6) 


or 


xdiere 


V  *  V  +  V 

p/c  p/B  pTC 


pTC 


X  r_  +  V. 


p/B 


B/C 


This  term  is  called  the  transport  velocity.  The  interpretation  of 
transport  velocity  defined  in  this  equation  is  such  that  V  is  still  the 
velocity  of  p  relative  to  B  and  is  the  velocity  in  C  that  p  would  have. 
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if  p  were  fixed  in  B.  Note  this  is  just  the  sum  of  the  translation  an( 
rotation  of  frame  B  relative  to  frame  C  if  the  point  p  is  considered  fixed. 

2.7. 4.2  Special  Acceleration.  By  taking  the  derivative  of  the  velocity,  as 
in  Equation  2.5,  and  applying  the  distributive  law  to  the  cross  product, 

c  c  c  c  c 

“  3t  'p/C  **  3t  'p/B  ’*clt  *^/C  ^  ^p/B  %/C^  3t  S/B  *  3t  '^B/C 

Now,  substituting  with  the  notation  for  acceleration  ^diere  possible, 

L  -  _  Cd  - 

^/c  =  at  '^P/B  +  ‘^/c  X  V3  +  0^/,  X  ^  rp/3  + 

The  two  remaining  terms  with  derivative  notation  should  be  recognized  as  appli¬ 
cations  of  Equation  2.4.  So,  substituting 


i 


\/c  =  +  %/C  X  V  ^3)  +  X  rp/3  +  X  (V /3  +  0^/^  X  +  A3 


/C 


Expanding  and  noting 


'^P/B  ^p/B 


-  V,  +  %/c  X  +  0^^^  X  X  V  X  X  )  + 

%/C 

Rearranging  and  combining  the  two  like  terms 

%/c  =  Vb  +  Vc  +  Vc  X  +  2Vc  X  V /3  +  X  X  r/3)  (2.7) 

Of  the  five  terms  remaining  in  the  acceleration  equation,  the  last  two 
have  descriptive  names. 


2“b/c  X  Vp^g  is  called  the  Coriolis  acceleration,  and 

%/c  X  %/c  X  Rp/B  called  the  Centripetal  acceleration. 

Ihe  terms  in  Equation  2.7  that  are  independent  of  the  motion  of  p 
relative  to  frame  B  are  called  the  transport  acceleration.  These  terms 
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provide  the  acceleration  in  frame  C  of  a  point  that  is  fixed  at  p  at  the 
instant  in  question.  Notationally,  the  transport  acceleration  is 


JL  —  ~  ^ 

^TC  “  \/C  ^/C  ^  ^P/B  %/C  ^  ^%/C  ^  ^p/B  ^ 


These  concepts  are  difficult  to  realize  until  a  few  problems  are 
attenpted. 

2. 7. 4. 3  Exaaple  Two  Reference  System  Problem.  The  angular  velocity  of  the 


FIGURE  2.16.  TWO  REFERENCE  SYSTEM  PROBLEM 


arm  ap  relative  to  the  disk  in  Figure  2.16  is  10  rad/sec,  shown  vectorally  in 
the  diagram  as  13^ ,  >diile  the  angular  velocity  of  the  disk  relative  to  the 
ground  is  5  rad/sec,  shown  vectorally  as  .  The  angular  accelerations  are 
zero.  Reference  B  is  attached  to  the  platform,  vdiile  frame  C  is  fixed  to  the 
ground,  three  feet  below  the  disk.  At  the  instant  in  question,  the  arm  ap  is 
in  the  vertical  position,  and  the  reference  axes  directions  coincide,  although 
displaced. 


2.31 


Find  the  velocity  and  acceleration  of  point  p  relative  to  the 
reference  frame  C. 

Using  Equation  2.6 

'p/C  “  ^p/C  ~  'p/B  ^/C  ^  S/B  ^B/C  (2.6) 

we  know  the  last  term/  V  .  =  0/  since  the  B  frame  is  only  rotating  relative 

to  C. 

=  5k  rad/sec  and  r^^^  =  3k  feet,  by  observation 
This  leaves  which  involves  angular  velocity  =  -lOi,  relative  to  B. 


\/B  “  “i  ^  '^p/B  “  (-lOi  X  3k)  =  (-30) (-j)  =  30j  ft/sec 

Substituting  all  the  parts  into  Equation  2.6 


=  30j  +  5k  X  3k  +  0  =  30j  +  15(k  X  k)  = 
For  the  acceleration,  the  general  expression  is  Equation  2.7 


30j  ft/sec 


\/c  '  \/B  +  Kc  +  ‘^/c  X  +  2(^/,  X  Vp/3  +  X  X  rp/3) 

The  only  unknown  terms  are  =  0  and  The  latter  is  a  centripetal 

acceleration  due  to  the  rotation  of  the  arm.  The  centripetal  acceleration  may 
be  arrived  at  in  several  different  ways, 

^  d  "  ■—  — ■ 

\/B  “  at  '^P/B  “  at  ^“1  ^  S/B )  =  «1  X  rp/3  +  X  rp/3 

=  0  +  X  (Wj  X  F/g)  =  (-lOi)  X  (30j) 

=  -300k  ft/sec^ 

Substituting  this  value  and  the  others  already  calculated 

=  -300k  +  0  +  0  X  3k  +  2(5k  X  30j)  +  5k  X  (5k  X  3k)  =  -300l  -  300k 


While  working  problems  v^ere  there  is  a  choice  of  axes,  be  careful  to 
choose  so  that  as  many  parameters  as  possible  are  equal  to  zero,  and  most 
inpartantly  so  that  the  axes  are  aligned  at  the  instant  in  question.  Also, 
\diether  a  reference  system  is  fixed  in  a  bocty  or  not  will  have  profound 
effects  on  the  velocities  as  seen  from  that  origin.  Try  to  place  yourself  at 
the  origin  of  a  system  and  visualize  the  velocity  and  acceleration  seen  to 
help  avoid  confusion.  Also  check  your  answers  to  see  if  they  are  logical, 
both  in  magnitude  and  direction.  The  right-hand  rule  is  essential. 

When  working  with  large  systems,  with  many  variables  it  becomes  necessary 
to  develop  a  shorthand  method  of  writing  systems  of  equations.  The 
development  of  matrix  algebra  is  the  solution. 

2.8  MATRICES 

An  m  X  n  matrix  is  a  rectangular  array  of  quantities  arranged  in  m  rows 
and  n  columns.  When  there  is  no  possibility  of  confusion,  matrices  are  often 
represented  by  single  capital  letters.  More  commonly,  however,  they  are 
represented  displaying  the  quantities  between  brackets;  thus. 


Note  that  a^^.  refers  to  the  element  in  the  ith  row  and  jth  column  of  [A]. 
Thus,  element  in  the  second  row  and  third  column.  Matrices  having 
only  one  column  (or  one  row)  are  called  col\jmn  (or  row)  vectors.  The  matrix 
[X]  below  is  a  column  vector,  and  the  matrix  [Y]  is  a  row  vector. 
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A  matrix,  ijnlike  the  determinant,  is  not  assigned  any  "value";  it  is 
simply  an  array  of  quantities.  Matrices  may  be  considered  as  single  algebraic 
entities  and  combined  (added,  subtracted,  multiplied)  in  a  manner  similar  to 
the  combination  of  ordinary  numbers.  It  is  necessary,  however,  to  observe 
specialized  algebraic  rules  for  combining  matrices.  These  njles  are  somewhat 
more  conplicated  than  for.  "ordinary"  algebra.  The  effort  required  to  learn 
the  rules  of  matrix  algebra  is  well  justified,  however,  the  simplification 
and  organization  vMch  matrices  bring  to  problems  in  linear  algebra. 


2.8.1 


TWO  matrices  [A]  =  [a.^]  and  [B]  =  [b.^]  are  equal  if  and  only  if  they 
are  identical;  i.e.,  if  and  only  if  they  contain  the  same  number  of  rows  and 
the  same  nianber  of  columns,  and  a.^  “  ^ij  values  of  i  and  j.  Thus, 


the  statement 


0  19 


is  ec[ui valent  to  the  statements 


=  2 
®12  =4 


2.8.2  Matrix  Addition 

Two  matrices  having  the  same  number  of  rows  and  the  same  number  of 
columns  are  defined  as  being  conformable  for  addition  and  may  be  added  by 
adding  corresponding  elements;  i.e.. 
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A  scalar  is  a  single  number.  A  matrix  of  any  shape  may  be  multiplied  Iw 
a  scalar  by  multiplying  each  element  of  the  matrix  by  the  scalar.  Ihat  is! 


2.8.4  Matrix  Multiplication 

Matrix  multiplication  can  be  defined  for  any  two  matrices  when  the  number 
of  coluTOS  of  the  first  is  equal  to  the  number  of  rows  of  the  second  matrix. 
This  can  be  stated  mathematically  as; 
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v^ere 


[A]  [B]  =  [C] 

ixm  inxj  ixj 


c.  . 
13 


Multiplication  is  not  defined  for  other  matrices. 

Equation  2.8  demonstrates  the  product  of  two,  2x2  matrices. 

[A]  [BJ  =  [C] 


2x2 

2x2 

2x2 

^11 

^12 

^12" 

S2' 

^21 

^2  2 

^21 

b2  2_ 

Si 

^2  2 

or  using  the  definition  of  multiplication, 


^11 

^12 

\i 

^12' 

^11 

+ 

^12 

i^21 

^11 

bl2 

+ 

^12 

^22 

^21 

^2  2 

>1 

^2_ 

^21 

bai 

+ 

S2 

^21 

^21 

^2 

+ 

S2 

^22 

This  situation  is  sufficiently  general  to  point  the  way  to  an  orderly 
multiplication  process  for  matrices  of  any  order. 

In  the  indicated  product, 


[A]  [B]  =  [c] 


the  left-hand  factor  may  be  treated  as  a  bundle  of  row— vectors. 


[A] 


®12  ^ 
•®2  2  ^ 
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and  the  rig^t-hand  factor  as  a  bundle  of  column  vectors. 


[B]  = 


Then, 


^  1  ^2  2^ 


and  ty  comparison  with  Equation  2.8 


^11  S  2 


Si  S2 


11  S  2  ^  ^1 1  ^1 


tSl  ^22  I  ^21  ^ 


v^ere,  by  definition, 


®i2^  Pill  f®ii  ®i2^  r^i 


®22^  I  111  ^Sl  ^22^  r^l 


(2.9) 


(2.10) 


[a^i  aj^j]  b^^  +  b^J 


ISi  S2I  ^2!  =  tSi  S2  ■»•  S2  ^22^ 


A  comparison  of  Equations  2.9  and  2.10  shows  that  if  the  rows  of  [A]  and 
the  columns  of  [B]  are  treated  as  vectors,  then  C.^  in  the  product  [C]  =  [A] 
[B]  is  the  dot  product  of  the  ith  row  of  [A]  and  the  jth  column  of  [B].  This 
rule  holds  for  matrices  of  any  size,  i.e.. 


2.37 


Matrix  multiplication  is  therefore  a  "row^on-column"  process: 


jth  colxjmn  ijth  element 


The  indicated  product  [A]  [B]  can  be  carried  out  only  if  [A]  and  [B]  are 
conformable;  that  is,  for  conformability  in  multiplication,  the  number  of 
columns  in  [A]  must  equal  the  number  of  rows  in  [B].  For  exanple,  the 
.expression 
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^11 

^12 

bx2' 

^21 

^22 

_ 

^21 

^22 

^32 

is  meaningless  (as  an  attempt  to  carry  out  the  multiplication  will  show) 
because  the  number  of  columns  in  [A]  is  two  and  the  number  of  rows  in  [B]  is 
three.  A  convenient  rule  is  this:  if  [A]  is  an  m  x  n  matrix  (m  rows,  n 
columns)  and  [B]  is  an  n  x  p  matrix,  then  [C]  =  [A]  [B]  is  an  m  x  p  matrix. 
That  is. 


(A] 


inxn 


tB] 


nxp 


[C] 


nxp 


Matrix  algebra  differs  significantly  from  "ordinary"  algebra  in  that 
multiplication  is  not  commutative.  In  general,  that  is. 


For  example,  if 


then 


[A]  [B]  t  [B]  [A] 


[B]  [A] 


2  -5 
4  2 
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Because  niultiplication  is  non-conimutative, 
the  product 


care  must  be  teOcen 


in  describing 


IC]  =  [A]  [B] 

to  say  that  [A]  "premultiplies"  fB],  or,  equivalently,  that  [B]  "post- 
multiplies"  [A]. 


2.8.5  The  Identity  Matrix 

The  identity  (or  matrix  [i]  occupies  the 

algebra  that  the  number  one  does  in  ordinary  algebra 
[A], 


same  position  in  matrix 
That  is,  for  any  matrix 


tl]  [A]  =  [A]  [I]  s  [A] 


Ihe  identity  (i] 

)  diagonal  and 


is  a  square  matrix  consisting  of  ones  on  the 
zeros  everyvdiere  else;  i.e.. 


principal 


1 


[I] 


0 

0 

* 

# 

0 


0  0  ...  0 
1  0  ...  0 

0  1  ...  0 

0  0  .  .  .  1_^ 


The  order  (the  number  of  rows  and  columns)  of  an  identity  matrix  depends 

entirely  on  the  requirement  of  conformability  with  adjacent  matrices.  For 
exanple,  if 
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2 


then 


[A] 


1  3 

0  1 


tU  [A] 


2  13 

1  0  1 


lA]  II] 


2  1 
1  0 


3 

1 


1  0 
0  1 
0  0 


0 

0 

1 


2  13 

10  1 


Thus,  the  "left"  identity  for  [A]  is  2  x  2  and  the  "right"  identity  for 
lA]  is  3  X  3;  however,  they  both  leave  [A]  unaltered. 


2.8.6  The  Transposed  Matrix 

The  transpose  of  [A],  labeled  [A]’’,  is  formed  by  interchanging  the  rows 
and  columns  of  [A].  That  is. 


1  2 

•  •  •  ®ln 

r  fj. 

^21 

♦  •  • 

®2 1  2 

*  •  *  ®2n 

^12  ®2  2 

•  •  •  a-2 

•  •  •  •  •  « 

•  •  •  •  •  • 

•  •  •  •  •  • 

1  0  •  •  •  3 

ml  in2  mn 

\ 

•  •  •  • 

•  •  •  • 

•  •  •  • 

•  •  •  cl 

mn 

The  jth  row  vector  becomes  the  jth  column  vector,  and  vice  versa.  For 
exanple. 
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3  1 

T 

3  2  0 

2  4 

14  2 

^0  2_ 

[2  5+j4 

T 

2  -10 

-10  3 

5+j4  3 

2.8.7  The  Inverse  Matrix 

Matrix  multiplication  has  been  defined;  it  is  natural  to  inquire  next  if 
there  is  some  way  to  divide  matrices.  There  is  not,  properly  speaking,  a 
division  operation  in  matrix  algebra;  however,  an  equivalent  result  is 
obtained  through  the  use  of  the  inverse  matrix. 

In  ordinary  algebra,  every  number  a  (except  zero)  has  a  multiplicative 
inverse,  a  ^  defined  as  follows:  A  quantity  a~^  is  the  inverse  of  a  if 


a  •  a  ^  =  a~^  a  =  1 


In  the  same  way,  the  matrix  [A]  ^  is  called  the  inverse  matrix  of  [A]  if 
[A]  [A]"^  =  [AJ"^  [A]  =  (I] 

The  symbol  1/a  is  normally  used  to  signify  a~^ .  Since  ordinary 
multiplication  is  comnaitative. 


(1/a)  .  (b)  =  (b)  •  (1/a)  =  b  -r  a 

for  any  mmiber  b.  The  use  of  the  division  symbol  (^)  in  this  instance  is 
useful  and  tinambiguous.  In  matrix  algebra,  however,  multiplication  is  not 
.commutative.  Therefore, 
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[A]"^  [B]  i:  [B]  [A]-^ 


and  the  expression 


[B]  ^  [A] 


cannot  be  used  since  it  nay  have  either  of  the  (unequal)  neanings  in  the 
previous  equation.  Instead  of  saying  "divide  [B]  by  [A],"  one  must  say  either 
"premultiply  [B]  by  [A]'^"  or  "postmultiply  [B]  by  [A]"\"  ihe  results,  in 
general,  are  different. 

2.8.8  Singular  Matrices 

Matrices  vdiich  cannot  be  inverted  are  called  singular.  For  inversion  to 
be  possible,  a  matrix  must  possess  a  determinant  not  equal  to  zero.  For 
example,  the  matrix 


2  1 

3  0 

4  5 

is  singular  because  it  is  not  square,  and  a  determinant  cannot  be  computed. 

The  matrix 

2  1 

4  2 


is  singular  because  its  determinant  vanishes. 

Matrices  \diich  do  possess  an  inverse  are  called  nonsingular. 
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2.9  SOLUTigy?  OF  LINEAR  SYSTEMS 

Consider  the  set  of  equations 


That  is. 


^21  +  ^22  *2  + 


*1  +  ®n2  *2  +  •  • 


[A]  [X]  =  [Y] 


+  a,  X 

In  n 


=  yi 

+  a2n  \  =  y2 


+  ®nn  K 
nn  n 


*  yn 


(2.11) 


Assuming  that  the  inverse  of  (A]  has  been  confuted,  both  sides  of  this 
equation  may  be  premultiplied  by  [A]"^ ,  giving 


[A]-'  lA]  [X]  =  (A]-"  [Y] 


-1 


From  the  definition  of  the  inverse  matrix. 


[1]  IX)  =  [AJ"^  [Y] 


frcan  which,  finedly. 


tx]  =  [A]-"  [Y] 


Thus,  the  system  of  Equation  2.11  may  be  solved  for  x^ ,  x^ ,  .  .  .  x^  by 
confuting  the  inverse  of  [A]. 


2.9.1  Conpiting  the  Inverse 

There  is  a  straightforward  four  step  method  for  confuting  the  inverse  of 
a  given  matrix  [A]: 


Step  1.  Conpite  the  determinant  of  [A].  This  determinant  is  written  as 
|Aj.  If  the  determinant  is  zero  or  does  not  exist,  the  matrix 
(A)  is  defined  as  singular  and  an  inverse  cannot  be  found. 
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step  2.  Transpose  matrix  [A].  oSie  resulting  matrix  is  written  [A]’^. 

Step  3.  Replace  each  element  a.^.  of  the  transposed  matrix  by  its  co¬ 

factor  A^  ^ .  This  resxilting  matrix  is  defined  as  the  adjoint  of 
matrix  [A]  and  is  written:  Adj  [A]. 

Step  4.  Divide  the  adjoint  matrix  by  the  scalar  value  of  the  determi¬ 
nant  of  [A]  vdiich  was  confuted  in  Step  1.  The  resulting  matrix 
is  the  inverse  and  is  written:  [A]~^ . 


This  procedure  can  be  summarized  as  follows:  To  calculate  the  inverse  of 
[A]  calculate  the  adjoint  of  [A]  and  divide  the  determinant  of  [A]  or 


[A]“^  =  tA] 

A. 

Example:  Find  [A]"^ ,  if 


[A]  = 


3 

1 

0 


2  0 

5  1 

2  -1 


Step  1.  Compute  the  determinant  of  [A].  Expanding  about  the  first  row 


A 

A 

A 

The  determinant  has  the  value 
Step  2.  Transpose  [A] 


3  2  0 

=  15  1 

0  2  -1_ 

=  3(-5  -2)  -2(-l  +  0)  +  0(2  -  0) 

=  -21  +2+0  =  -19 

—19;  therefore  an  inverse  can  be  computed. 


[A]^  = 


3 

2 

0 


1  0 

5  2 

1  -1 
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step  3:  Replace  each  element  a^ .  of  [A]’^  by  its  cofactor  A.  .  to  deter¬ 
mine  the  adjoint  matrix^  Note  that  signs  alternate^ trom  a  posi¬ 
tive  Ajj 


Adj  [A]  = 


2 

-1 

0 

-1 

0 

2 


2 

-1 

0 

-1 

0 

2 


-7  2  2 

1  -3  -3 

2  -6  13 


Step  4:  Divide  by  the  scalar  value  of  the  determinant  of  [A]  \diich  was 
confuted  as  -19  in  Step  1. 


[A]-' 


-7  2 

1  -3 

2  -6 


2.9.2  Product  Check 

From  the  definition  of  the  inverse  matrix 


[A]-"  [A]  =  [I] 


This  fact  may  be  used  to  check  a  coitpated  inverse.  In  the  case  just  completed 


2.46 


Since  the  product  does  come  out  to  be  the  identity  matrix,  the 
conputation  was  correct. 

2.9.3  Exanple  Linear  System  Solution 

Given  the  following  set  of  simultaneous  equations,  solve  for  ,  x^ ,  and 

X3. 

3Xj  +  2X3  -  2X3  =  y^ 

-X3  +  Xj  +  4X3  =  y^ 

2X3  -  3X3  +  4X3  =  y3 

This  set  of  equations  can  be  written  as 

[A]  [X]  =  [Y] 

[X]  =  [A]-'  [Y] 


Thus,  the  system  of  Equations  2.12  can  be  solved  for  the  values  of  x^ , 
Xj ,  and  X3  hy  computing  the  inverse  of  (A] . 
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lA]  [X] 


Step  1: 


Step  2: 


Step  3: 


=  [Y] 


3  2-2 

-114 
2-3  4 


Yi 

^2 

ya 


CcMipite  the  determinant  of  [A].  Expanding  about  the  first  row 

A  =  3  (4  +  12)  -2  (-4  -8)  -2  (3  -  2) 

A  =  48  +  24  -  2  =  70 


Transpose  (AJ 


lAr 


3 

2 

-2 


-1  2 

1  -3 

4  4 


Determine  the  adjoint  matrix  hy  replacing  each  element  in  [A]^  by  its 
cofactor 


Adj  lA]  = 


1  -3 

2  -3 

2  1 

4  4 

-2  4 

-2  4 

Adj  [A]  = 


-1 

4 

-1 

1 

16 

12 

1 


3 

-2 


|2 

-3| 

-2 

10 

16 

-10 

13 

5_ 

3 

-2 


-1 

4 


3  -1 

2  1 
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step  4:  Divide  by  the  scalar  value  of  the  determinant  of  [A]  which  was 
conputed  as  70  in  Step  1. 

16  -2  10 

[A]"^  =  ^  12  16  -10 

1  13  5 

Product  Check 


[A]-'  [A]  =  [IJ 


Since  the  product  in  Equation  2.13  is  the  identity  matrix,  the 
computation  is  correct.  The  values  of  ,  x^ ,  and  x^  can  now  be  found  for  any 
y^,  Yj,  and  by  premultiplying  [Y]  by  [A]“\ 
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8 


For  exanple,  if  ■=  1,  Yj  =  13,  and  = 


16  -2  10 

1 

^2 

II 

12  16  -10 

13 

*3 

1  13  5 

8 

--  — 

Xj  ®  (16  -  26  +  80)  =  =1 

Xj  =  ^  (12  +  208  -  80)  =  ^  =  2 

X3  =  ^  (1  +  169  +  40)  =  Ir  =  3 

Solution  of  sets  of  simultaneous  equations  using  matrix  algebra 
tecdiniques  has  wide  application  in  a  variety  of  engineering  problems. 
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PROBLEMS 


2.1 


Is 


-  V?  -  ^/T-  -vT- 

V  -  i  +’-J  j  k 


a  unit  vector? 


2.2  Find  a  unit  vector  in  the  direction  of 

A  =  2I  +  3  j  -  k 

2.3  Are  the  following  two  vectors  equal? 

A  =  2i  +  3j  -  k 
i  =  4i  +  6 j  -  2k 

2.4  The  following  forces  measured  in  pounds  act  on  a  body 

=  2I  +  3j  -  5k 
Fj  =  -5i  +  j  +  3k 
F3  =  i  -  2J  +  4k 
F^  =  4i  -  3 j"  -  2k 

Find  the  resultant  force  vector  and  the  magnitude  of  the  resultant  force 
vector. 

2.5  If  A  =  3i  -  j  -  4k 

B  =  -2i  +  4 J  -  3k 
C  =  i  +  2j  -  k 

Find 

2A  -  B  +  3C 
|a  +  B  +  cj 
I  3A  -  2B  +  4C  I 
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2.6  Hie  position  vectors  of  points  P  and  Q  are  given 

r^  =  2i  +  3j  -  k 
r^  =  4i  -  3j  +  2k 

Determine  the  vector  from  P  to  Q  (PQ)  and  find  its  magnitude. 


2.7  Find  A  *  B  using  A  and  B  from  Problem  2.5. 


2.8  Given  A  =  2l  +  3]  -  Tc 

S  *  41  +  6J  -  2Tc 

a.  Find  A  •  B  _  ^ 

b.  Find  the  angle  between  A  and  B. 


2.9  Evaluate 

j.(2l  -  3j  +  k)  = 
(2T«J)»(3i  +  Ic)  = 


2.10  If  A  =  3l  -  J  -  4k 

B  =  -2i  +  45  -  3lc 

Find  A  X  B 

2.11  Determine  the  value  of  "a"  so  that  A  and  B  below  are  perpendicular. 


A  =  2i  +  aj  +  it 

B  =  4l  -  2  J  -  21^ 

2.12  Determine  a  xmit  vector  perpendicular  to  the  plane  of  A  and  B  below. 

A  =  2i  -  6T  -  3lc 

B  =  4i  +  3]  -  1c 

2.13  If 

A  =  2I  -  35  -  ic 

B  =  T  +  4j  -  2lc 

Find 


A  X  B 
B  X  A 


2.52 


«nd  _  _  _  _ 

(A  +  B)  X  (A  -  B)  (the  quick  way  ijsing  vector 
algebra) . 

2.14  Evaluate 

a.  2j  X  (3i  -  4k) 

b.  (I  +  2j)Xk 

2.15  The  aircraft  shown  below  is  flying  around  the  flagpole  in  a  steady  state 
turn  at  a  true  velocity  of  600  ft/sec.  The  turn  radius  is  6,000  ft. 
What  is  turn  rate  w  expressed  in  unit  vectors  (i,  j,  k)  of  the  XTZ  system 
shewn? 


2.16  For  the  same  aircraft  and  conditions  as  Problem  2.15,  vhat  is  turn  rate 
expressed  in  unit  vectors  (i,  j,  k)  of  the  xyz  system  shown? 

2.17  Given 

r  “  t^i  -  6tj  +  6k 

Find  r  with  respect  to  the  axis  system  xyz  ^ich  has  i,  j,  k  as  its  unit 
vectors.  Is  r  a  velocity? 
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2.18  If  the  3^  system  in  Problem  2.17  is  rotating  at  3i  +  2j  -  k  rad/sec  wifll 
respect  to  another  system  XYZ,  find  T  with  respect  to  XYZ.  Is  r  the 
velocity  of  the  point  \diose  radius  vector  is  r  with  respect  to  XIZ?  What 
system  is  the  answer  of  this  problem  referred  to? 

2.19  A  flyvdieel  starts  from  rest  and  accelerates  counterclockwise  at  a 
constant  3  rad/sec^ .  After  six  seconds  the  point  P  on  the  rim  of  the 
%dieel  has  reached  the  position  shown  in  the  sketch.  What  is  the  velocity 
of  point  P  with  respect  to  the  fixed  XYZ  system  shown? 


PAPER 


Y 


2.20  If 


A  «  3t^i  -  tj 
B  -  -6ti  +  tk 


Find 

vectors. 


d(A»B)/dt  relative  to 
Is  the  answer  a  vector? 


the 


system  having  i. 


and  k  as  its  unit 
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21  A  small  body  of  mass  m  slides  on  a  rod  which  is  a  chord  of  a  circular 
wheel  as  shown  below.  ihe  wheel  rotates  about  its  center  with  a 
clockwise  velocity  4  rad/sec  and  a  clockwise  angular  acceleration  of  5 
rad/sec^ .  Hie  body  m  has  a  constant  velocity  on  the  rod  of  6  ft/sec  to 
the  right.  Relative  to  the  fixed  axis  system  XY  shown  below,  find  the 
absolute  velocity  and  acceleration  of  m  when  at  the  position  shown. 
Hint:  Let  xy  system  rotate  with  the  disk  as  shown. 


(0  -  -4k 
w  -  -5k 


! 

I 
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2.22  A  small  boy  holding  an  ice  cream  cone  in  his  left  hand  is  standing  on  thU 
edge  of  a  carousel.  The  carousel  is  rotating  at  1  ra^sec  counterclock¬ 
wise.  As  the  boy  starts  walking  toward  the  center  of  the  \dieel,  \diat  is 
the  velocity  and  acceleration  vector  of  the  ice  cream  cone  relative  to 
the  ground  Xi7 

Hint:  Let  xy  be  attached  to  the  edge  of  the  carousel. 


Bey's  velocity 
Bey's  acceleration 
Carousel's  acceleration 


2  ft/sec  toward  center 
1  ft/sec^  toward  center 
1  rad/sec^  counterclockwise. 


2.23  Solve  the  following  equations  for  ,  and  by  use  of  the  inverse 

matrix. 


X,  +  X3 


2X2  +  2X3  »  1 


-Xj  +  Xj  +  2X3  «  3 
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2.24  For  a  -  i  let 


Conpite 

a. 

b. 

c. 

d. 

e. 

f. 

g- 

h. 

i. 
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[A]  = 


[C]  = 


2 

3 


2  -1 

-3  2 


[X]  = 


IC]  -  lA] 

[A]  [B] 

([A]  [B])  [Y] 
[A]  ([B]  [Y]) 
[A]  [C] 

IC]  [A] 

[Xf  [B] 
txf  ([A]  [X]) 
IX]  iDf 


If  lA]  = 


2.26  Find  x,  y,  and  z 


-1 

1 

0 


IB]  = 


ID]  = 


lY]  = 


4 

-1 

4 

2 

Yi 

Ya 

Ys 


-4x  -  3y  -  3z  =  1 
X  +  z  =  1 

4x  +  4y  +  3z  =  1 


3 

1 


Find  IA]IA] 
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2.27  Read  the  question  and  circle  the  correct  answer.  True  (T)  or  False  (F); 


T  F  A  vector  is  a  quantity  whose  direction  and  sense  are  fixed,  but 

vdiose  magnitude  is  unspecified. 

T  F  A  scalar  is  a  quantity  with  magnitude  only. 

T  F  The  magnitude  of  a  unit  vector  is  one. 

T  F  Zero  vectors  have  any  direction  necessary. 

T  F  A  free  vector  can  be  moved  along  its  line  of  action,  but  not 

parallel  to  itself. 

T  F  Free  vectors  may  be  rotated  without  change. 

T  F  A  3  X  2  matrix  can  pre-multiply  a  2  x  4  matrix  and  the  result 

will  be  a  3  X  4  matrix. 


T  F  A  3  X  2  matrix  can  post-multiply  a  2  x  4  matrix  and  the  result 
will  be  a  3  X  4  matrix. 


T  F 

T  F 


Multiplying  a  matrix  by  a  scalar  is  the  same  as  multiplying  its 
determinant  by  the  same  scalar. 

I 

Identity  matrices  are  always  square.  " 


T  F 


0 

1 


T  F  Both  matrices  in  the  preceding  question  are  identity  matrices. 
T  F  Singular  matrices  can  be  inverted. 


T  F 
T  F 
T  F 

T  F 
T  F 


The  determinant  of  a  non-singular  matrix  is  zero. 
Inverting  a  matrix  is  a  straightforward  process. 


A  X  B 


i  j  k 

Ai  ^  ^3 


B. 


B, 


B, 


1  2  “3 

The  determinant  of  any  matrix  can  be  calculated. 

The  value  of  a  determinant  depends  upon  \diich  row  or  coltmin  it 
was  expanded  about. 
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T  F  Velocity  is  the  time  rate  of  change  of  a  velocity  vector. 

T  F  Acceleration  is  the  time  rate  of  change  of  a  velocity  vector. 

T  F  Acceleration  has  to  be  expressed  in  (referred  to)  unit  vectors 

of  an  inertial  reference  system. 

T  F  Bodies  moving  with  pure  translation  only  do  not  rotate. 

T  F  Reference  systems  are  considered  to  be  non-deformable  rigid 

bodies. 

T  F  [A]  [B]^  =  [B]  [A],  if  the  two  matrices  are  conformable  for 

multiplication  on  the  left  hand  side  of  the  equation. 

T  F  Ivl  =  |-v| 

T  F  |S|  »  I  i 

T  F  The  magnitude  of  vIaI  is  equal  to  B/1b| 

T  F  2(3A)  =  5A 

T  F  if  j/  and  k  are  orthogonal. 

T  F  IpqI  is  the  distance  between  points  P  and  Q. 

T  F  A*B=B*A 

T  F  If  A  •  B  is  zero  and  neither  A  nor  B  are  zero,  then  A  and  "i 
must  be  parallel. 

T  F  i  -7  =  1 

T  F  AXB=BXA 

T  F  AXB^A^Bj+AjBj+AjBj 

2.28  Define: 

Determinant 

Vector 

Scalar 

Free  vector 

Bound  vector 

Velocity  vector  of  a  particle 


unit  vector 


Zero  vector 
Parallel  vectors 
PositicHi  vector 
Matrix 

Square  matrix 
Colunn  vector 
Bow  vector 
Matrix  equality 
Matrix  conformability 
Matrix  non-commutativity 
Identity  matrix 
Transposed  matrix 
Singxilar  matrix 

2.29  Find  the  VW's  velocity  and  acceleration  vectors 


2.30  Find  a  unit  vector  parallel  to 

A  «  2i  -  3 j  +  6k 


2.31  What  is  the  magnitxide  of  the  following  vector? 


A  «  2i  +  3j  +  6k 


2.32  Is  the  axis  shown  a  "right-handed"  axis  system? 


2.34  Add  the  following  vectors 

A  =  3i 
B  «  4k 
C  -  1/2 j 
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2.39  The  angular  velocity  of  a  rotating  rigid  body  about  an  axis  of  rotation 
is  given  by  w  ®  4i  +  2j  +  k.  Find  the  linear  velocity  of  the  Point  P  on 
the  body  wtose  position  vector  relative  to  a  point  on  the  axis  of 
rotation  is  T  -  2j  +  2k. 
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2.40  The  T-38  shown  is  in  a  right  continuous  roll  at  2  racVsec  %^ile  travelin' 
at  480  ft/sec.  Find  the  velocity  of  the  wingtip  light  with  respect  to 
the  axis  shown. 


Y 


XYAXIS  IS  INERTIAL 
(i.e.  FIXED) 


12.5 

FEET 


2.41  The  particle,  P,  is  following  a  path  described  ;  x  *  6t^ ,  y  *  t  +  1,  z 
«  t^ .  Find  the  velocity  and  acceleration  of  P,  with  respect  to  the  axis 
shown. 
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If  A  -  sT  -  j  +  2k 

and  B  «  -i  +  3j 

-k,  find 

a.  IaI 

f. 

A  X  B 

b.  |f|  _ 

9- 

Unit  vector,  a,  parallel  to  A 

C.  A  +  B 

1  1 

h. 

a  •  A 

d.  1  A  B  1 

i. 

A  — 

a  •  B 

e.  A  •  B 

2.43  Ihe  shaft  is  rotating  counterclockwise  around  the  cone  in  the  XZ  plane  at 
5  ract/sec  and  accelerating  at  3  rad/sec^ .  The  \dieel  is  rotating  as  shown 
at  200  rad/sec  and  decelerating  at  50  rad/sec^ .  Find  the  velocity  of 
point  P  with  respect  to  reference  system  C  at  the  instant  shown.  Hint: 
Let  X  be  fixed  in  the  shaft,  and  xz  and  XZ  planes  remain  coplanar. 


1.5FT 

RADIUS 
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2.44  If  [A]  = 

1 

-2 

3 

and 

IB]  = 

1 

-1 

2 

1 

0 

0 

3 

Find:  [A] 

_2 

1_ 

[B]  and 

[B]  [A] 

2.45  If  [A]  - 

1 

2 

and 

[B]  * 

** 

1 

0 

1 

3 

■P  m 

Find  [A]  [B]  and  [B]  [A] 
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2.47 


2.48 


3 

0 

1 

-1 


+(3) 


If 


Find 

2.49  If 

Find  k. 


X  1 

2  y  +  z 
X,  Y,  and  z. 


2 

-a 


=  k 


2.50  Ccxtpate  the  inverse  of 

2  -1 

-3  2 


-1 

2 

1 

2 


-4 

2a 


3 

0 

0 

1 


4  y  -  z 
2  5 
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2.51  Conpjte  the  inverse  of 

3  2  1 

15  4 

6  4  2 

2.52  Coipite  the  inverse  of 

2  4  1 

-10  2 

11  1 


2.53  For  v^t  value  of  y  is  this  matrix  singular? 

“T 

1  3  y 

2  0-1 

1  1  -y 


2.54  Find  the  determinant  of 

6  0  0  0 
8x00 


0  0 
0  0 
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2.55  If  [A] 


2 


4 


1 


-10  2 
111 

Find  [A]  [Aj 

2.56  If 

X  +  2y  +  32  «  dj 

4x  +  5y  +  62  »  dj 

7x  +  8y  +  92  *  83 

Find  X,  Yt  z  for  any  value  of  a^ ,  a^ ,  and  aj . 
Find  X,  y,  and  z,  vdien  a^  *  1,  a^  *  2,  and  83  «  3. 


2.57  The  mechanism  shown  below  vibrates  £dx}ut  its  ec[uilibrium  position,  E.  At 
the  instant  shown  block  A  has  a  velocity  of  5  ft/sec  to  the  right  and  is 
decelerating  at  4  ft/sec^  to  the  left.  The  bob  B  in  its  counterclockwise 
mortion  maintains  a  constant  angular  velocity  |  si  of  5  rad/sec.  Calcul^^^W 
the  velocity  and  acceleration  of  the  bob  relative  to  the  given  XY  systSi 
at  the  instant  shown.  Hint:  Let  the  xy  axis  be  fixed  to  the  block  A. 
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2.58  Capt.  Marvel,  US  Amy,  is  performing  a  loop  with  an  angular  velocity  jco] 
of  1  ra^sec  in  his  Ruey  Cobra  to  roll  in  on  a  target.  At  the  top  of  the 
loop,  the  leading  rotor  blade  is  just  parallel  with  the  helicopter's 
centerline.  The  rotaticm  of  the  rotor  is  |  is  3  rad/sec  counterclockwise 
as  viewed  from  the  top  of  the  helicopter.  At  this  instant,  \diat  is  the 
velocity  of  the  leading  rotor  blade  tip?  If  Capt.  Marvel  were  to  raise 
the  collective  and  accelerate  the  rotor  speed  3  racl/sec^,  this  would 
accelerate  an  angular  velocity  of  his  loop  1  raci/sec^ .  What  would  the 
acceleration  of  the  leading  rotor  blade  tip  be?  Hint:  Let  the  xyz 
system  be  attached  to  the  helicopter  rotor  path  plane  as  shown. 

radius  of  loop  -  1,000  ft 
radius  of  rotor  path  plane  ■  10  ft 


Hint:  Let  the  YZ  and  yz  planes  remain  coplanar. 
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ANSWERS 


2.1  Yes,  magnitude  of  V  =  1 

2  2  3j  -k 
(14) 

2.3  No,  B  =  2A 

2.4  Fy  »  21  -  J 

l?xl.  - 

2.5  111  -  8k;  ^^93  ;  «f398 

2.6  PQ  =  2l  -  6 j  +  3k 

IpqI  =  >J49 

2.7  A  •  B  =  2 

2.8  A  •  B  =  28;  4>  =  0 

2.9  -3;  undefined  i  •  j  =  0 

2.10  A  X  i  =  19i  +  17j  +  loic 

2.11  a  =  3 

3*  2  <  6t 

2.12  u  =  ■;^l  -  +  ^k 

2.13  A  X  B  =  lOr  +  3j  +  Ilk;  B  X  A  =  -loT  - 
(A_+  B)JX.  (A  -  B)  =  -20r  -  6j  -  22k 

2.14  -8i  -  6k;  +2i  -  j 

2.15  “  ®  j 


i  -  Ilk; 
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2.16  (0  =  .07j  -  .07k 

2.17  r  =  3t^r-6j 

2.18  (r)^y^  =  (3t^  -  6t  +  12)i 

-(t^  +  24 )j 
-(2t^  +  18t)lc 
in  XYZ  system 


2.19  V 

xy  z 

=  -108j 

ft/sec 

2.20  d(A  •  B)  =  -54t^ 

— at — 

2.21 

P/C 

=  181  + 

6j; 

^/C 

39i  -  88. 5 j 

2.22  V,,. 

P/C 

=  -lOi  - 

•  j; 

^/C  ” 

-loj  -  7I 

2.23  = 

=  -1/4;  > 

‘2 

II 

1 

X3  =  9/4 

2.25 

5  -3 

1 

2  1 

4 

3  -1 

2 

2.26  X  = 

-10;  y 

= 

2;  z  = 

11 

2.29  V  = 

16ti;  a 

= 

161 

2.30  ,.a  = 

-  2i  -  3r  + 

-  7I  7: 

^k 

2.31  |a| 

=  7 

2.33  a  = 

12k 
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2.34  A+B  +  C  =  3l  +  4k  +  l/2j 

2.35  A  +  B  =  8i  +  3j;  <|>  =  20" 


2.36  ♦  =  27.6" 


2.37  A  X  B  =  -28i 

2.38  A  X  B  *  -2i  -  5j  +  2k 

2.39  Vp  =  61  -  7 j  -  10k 

2.40  V  =  480l  +  25k 

2.41  V  =  12tr  +  J  +  3t^k;  a  =  12i  +  6tk 

2.42  a  f  -5i  +  j  +  8k _  _ 

b  sTH  g  (3//l4)i  -  (l//l4)j  +  (2A/l4)k 

c  2T  +  2j  +  k  h  ^/l4 

d  3  i  -8/yiI 

e  -8 


2.43  V, 


_P/C 

^/c 


2.44 


[A]  [B] 


2.45 


lA]  tB] 


2.46 


300 j  +  7.5i  -  15k 
-70.51  -  75 J  -  60046.5k 


7  -4 
2  1 


[B]  [A] 


L. 


tB]  [A] 


9 

2 


-1  -3 

=  6  3 

4  -3 

Cannot  do 


3 

0 

6 
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2.47 


4  3 

5  5 


2.48  x-4;y«3;z»2 


2.49  k  -  -1/2 


2.50 


(A] 


-1 


2  1 

3  2 


2.51  Ik>  Inverse 
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[A] 


-1 


1 

1 


2.53  y.  -  1/4 

2.54  72 

2.55  |”l 
0 
2 

2.56  No  Solution 


-2 

3 


9 

-2 

5 


-3 

1 


8 

-5 


11 

1 

4 


2.57  «=  13.66i  +  5j;  A^^^  =  -29i  +  43. 3j 


2.58  V, 


p/c 


-30i  -  1000 j  -  10k;  =  -30i  -  900 j  -  1010k 
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